First Order Transition in the Ginzburg-Landau Model 
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The d-dimensional complex Ginzburg-Landau (GL) model is solved according to a variational 
method by separating phase and amplitude. The GL transition becomes first order for high super- 
fluid density because of effects of phase fluctuations. We discuss its origin with various arguments 
showing that, in particular for d = 3, the validity of our approach lies precisely in the first order 
domain. 
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Ginzburg-Landau free energy functionals involving a 
ri-component space and time dependent field have been 
widely used in order to describe different types of phase 
transitions in a semi-phenomenological way. The case 
n = 2, corresponding to a complex field ip, applies 
in particular to superconductivity, superfluidity, metal- 
insulator transitions or to magnetic systems with mo- 
ments that are confined to a plane (XY moments). For 
such particular applications, the coefficients determining 
the functional can be derived from appropriate micro- 
scopic models. For superconductors, this was originally 
done by Gorkov Q| and has been refined since by numer- 
ous authors ||]. In its simplest form, the GL functional 
involves a time independent field and thus describes (clas- 
sical) thermodynamic and static phenomena of super- 
conductors. For the description of dynamic phenomena, 
charging effects 1^ or pairing fluctuations in strong cou- 
pling superconductors ||^ , one needs the generalisation to 
a time dependent field. 

An important aspect of the case n = 2 is the interplay 
between variations of amplitude and phase (j) of the 
corresponding complex field ■(/' = IV'le*"^. In various ap- 
proximate treatments, such as the mean field approxima- 
tion or the Hartree decoupling of terms involving higher 
powers of '0, the field is treated as a whole, without sep- 
arating amplitude and phase. More accurate studies of 
the static GL problem, like the renormalisation group ap- 
proach focusing in particular on the region near the 
phase transition, show that the amplitude has no crit- 
ical behaviour and is irrelevant at the transition. The 
phase transition scenario should then correspond to the 
one of the XY model with the same dimensionality. In 
the framework of the e-expansion IJ^, for d = 3, the 
transition is second order and seems to have the same 
critical exponents as the XY model. On the other hand, 
Bormann and Beck ||^ have shown that amplitude fluc- 
tuations, even though not being critical by themselves, 
might alter the cooperative phenomenon occurring with 
phases, in particular in dimension 2. Like the XY model, 
corresponding to a fixed value of the amplitude, the 2(iGL 
model can be mapped onto a Coulomb gas describing 
vortex-antivortex pairs. As soon as one allows for am- 
plitude variations, these topological excitations become 
energetically more favorable. Taking into account gaus- 
sian amplitude fluctuations, Bormann and Beck j|] have 
shown that the system may be driven into a regime where 



- according to Minnhagen's phase diagram |^ - a first or- 
der transition replaces the usual Kosterlitz-Thouless sce- 
nario. 

As far as superconductors are concerned, BCS theory 
pO| predicts that the transition between the normal state 
and the superconducting state is a second order phase 
transition. However, it is well known that fluctuations 
can change the order of the transition. For example, fluc- 
tuations of the magnetic field change the GL-BCS transi- 
tion to a first order transition for type I superconductors 
The three state Potts model in two dimensions is 
an opposite example: mean field theory predicts a first 
order transition, whereas the actual transition is contin- 
uous. So the question of the order of the transition in the 
GL model, as well as the more detailed mutual influence 
between phase and amplitude, is still an open problem. 

The aim of this letter is to show the reciprocal influ- 
ence between phase and amplitude by separating self- 
consistently, from the outset, the GL functional into two 
parts: the amplitude part and the phase part. 

According to Ginzburg-Landau theory, we define the 
effective hamiltonian functional 



(1) 



where a, h and 7 are coefficients independent of the tem- 
perature derived from a microscopic model, t = T/Tq — 1 
is the reduced temperature and Tg is the mean field crit- 
ical temperature. We now introduce the amplitude ji/'l 
and the phase (j> of the field tp = |'(/;|e"^. On the lattice, 
with lattice spacing s, we normalize the hamiltonian by 
setting i?^ = |V'P/(a/6), u — r/Co, where — 7/a is the 
mean field correlation length at zero temperature. The 
normalized hamiltonian is then: 



N 
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where /i := [1 " cos(0i - 0^ )] and 



N 



where Vi?i :— ^ ^j)- J points to the nearest 

neighbours of i and ej is a unit vector in direction j. We 
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have set R^Rj = Ri{{Rj - Ri) + Ri) « i?f in the XY 
part of the hamiltonian. Indeed, the term Ri{Rj — Ri) 
is not important for the present discussion because our 
approach will be reliable when amplitude fluctuations are 
small compared to phase fluctuations. Only two param- 
eters, that are in competition, remain: 



Vo := -j—Tl£ 

Kb 



d~2 



a controls amplitude fluctuations. Vq corresponds to the 
zero temperature phase stiffness. Vq is proportional to 
the supcrfluid density ajh and controls the general critical 
behaviour. When Vb/To is large, the critical region is 
small and the material behaves according to mean field 
theory (as BCS superconductors). When Vq/Tq is of the 
order of 1 , phase fluctuations become very large and give 
an upper bound for the critical temperature [ p^ . One 
has to bear in mind that Vq is not a constant that only 
normalizes the temperature T in the Boltzmann factor 
because of the occurrence of the temperature t = T/Tq — 1 
inside the hamiltonian 

Let us write the partition function in polar coordinates: 

Z = / D(i> \ DR e'l^"-" 
Jo Jo 

where the effective hamiltonian H^ff is 



T 



(3) 



We keep the factor R of the Jacobian, R Dcf) DR, and 
take it in the exponential giving a contribution logR to 
the potential. This incorporates the fact that small val- 
ues of R have a small statistical weight, due to the volume 
element in phase space, into the Boltzmann factor of the 
canonical ensemble. We compute now the partition func- 
tion by integrating only the phase. For this purpose, we 
will drop the integration on Ri and search for the mini- 
mum of the free energy with respect to R. The partition 
function becomes then: 



D(j) e 



The minimum of the free energy F — —ksT log Z is given 
through the equation 5F/6Ri — 0. Assuming that the 
gradient of the amplitude is zero, we have 



a (tR^ + 



ZVq 







(4) 



where we have multiplied the equation by R/2. f{K) = 
(fi) is the expectation value within the XY model of fi 
with a dimensionless coupling constant K ~ Y^R^. Al- 
though K has the same value at each lattice site and does 
not explicitely adapt to the vortex structure of the phase 
field, the average energy of the latter still determines the 
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FIG. 1. Reduced mean amplitude as a function of 
the reduced temperature for d = 3 for different values 

of d. The slanting dashed line shows the XY transition 
ending at — uJ^d = 2.2 (a oo). The vertical dashed 
lines indicate the temperature where the first order transition 
arises. The free energy for cr = 2 is also shown (in 

arbitrary units). 



value of K through the minimalisation of the free energy 
F. In this work, we take Monte-Carlo simulations to eval- 
uate the function f{K) which is just the energy of the 
XY model The critical temperature Tc is reached 

when the coupling K equals the critical constant 1/w: 



l = uK,^u^{R^){T,) 



(5) 



where a; = 0, 0.9, 2.2 (for d = 1,2,3 respectively) is the 
pure XY critical temperature on a square lattice with 
unit coupling constant K — \/T. 

The solution of equation (^) is plotted in figure 1 for 
d = 3. We use a normalization which is independent of 
the temperature t: with R^ = R^/{-t), f ^ T/{-tVo) 
and (T = —ta. The free energy F, calculated as: F ~ 
J dR dF/dR, has three branches. The crossing point of 
the two lower branches determines the location of the 
first order transition, marked by a vertical dashed line. 
The consequence is a first order transition for 



= -atr < 4.5, 



(6) 



For d = 2, 
< 1.25. For 



i.e. the mean amplitude R^ jumps at Tc 
we have also a first order transition for ct, 
d ~ I, there is no transition. 

A good approximation for the critical temperature due 
to phase fluctuations is given by (P) with (i?^) = —tc- 



(7) 



Combining this last equality and (g) , the first order tran- 
sition occurs then for (tTc/Vo < C where C = u) 4.5 « 9.9 
for d — i. 

What are the influence of amplitude fluctuations 
around the saddle point that we found? Do they de- 
stroy the first order transition or not? To answer to this 
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question, we include harmonic amplitude fluctuations in 
our computation. We separate phase and amplitude, in- 
troduce three variational parameters and derive two self 
consistent equations for these parameters. 

The effective potential Uefj ■= a [tRf + ^Rf) ~ 
^ log Ri of the effective hamiltonian ^ has only one 
minimum for all temperatures. So one can expect that 
a gaussian approximation for the amplitude gives a good 
approximation for all temperatures. We want also to get 
an hamiltonian with no direct coupling between phase 
and amplitude, but with effective constants that keep the 
memory of their interaction. The idea is then to separate 
phase and amplitude as: 



(8) 



Therefore, we set the trial hamiltonian: 

Ht[R, 4>]^Y.'^B [R, - Rof + i (Vi?,)' + A /,] (9) 

i 

Using the Bogoliubov inequality, we have F < Ft + 
{Heff — Ht)^ — F, where F is the free energy and (...)< 
is the canonical average with respect to Ht- The right 
hand side is thus to be minimized with respect to the 
constants A, Rq and B to give the best approximation 
of F. We introduce also the local amplitude fluctuation 
rji = Ri — Rq. The derivative of F with respect to these 
parameters gives three equations: 

A = (10) 
H + 3(ry2)) + (/)] i?2 + ^i?4 _ ^ (11) 

a{t + Siv')) + if) -B + iaRl ^ q (12) 

where all indices t and i are dropped, (logi?.;) is com- 
puted by a cumulant expansion. Equations ( pi] ) and ( [T^ ) 
are the central result of this work. The mean square am- 
plitude fluctuation is 

fc|<A 

where we extend the lower bound —Rq to — oo assuming 
that amplitude fluctuations are small. A is the reduced 
cut-off parameter and is computed on the first Brillouin 
zone. 

Equation (|ll|) reduces to equation (jj) when the addi- 
tional amplitude fluctuations are zero: (77^) — 0. The 
phase diagram is approximately the same as the one 
without amplitude integration, except that the critical 
temperature is smaller due to the additional harmonic 
amplitude fluctuations. The transition is first order for 
cr < 9.9 Vb/Tc for d = 3. 

In figure |^, we compare the simulations of Nguyen and 
Sudb0 jisl l with our phase-amplitude separation method 
for a corresponding a/ (Vq/To) = 18 (continuous regime). 
We see that the simulated (i?^) shows an inflexion point 
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FIG. 2. Expectation value of for d = 3 in the 
continuous regime. The thin line shows the simulation. 
We set the analytical critical temperature tc ~ —0.38 and the 
one from the simulations ts ~ —0.45. 



at Tc as the analytical (i?^) does. This is perhaps the 
forerunner for a jump at a smaller a. Both methods 
agree with the fact that there is no first order transition 
at large a. To our knowledge, accurate simulations for 
smaller a are not avalable. 

A remarkable thing is the fact that a negative cr pro- 
duces a first order transition ( ||l^ p. 340). Therefore, the 
critical line which is believed until now to separate the 
region a < and cr > is pushed in the region with cr > 
by phase fluctuations. The limiting case cr < is then 
consistent with a first order transition for a < 9.9 Vq/Tc- 
We mention here that a variational approximation due 
to Halperin et al |l^ for gauge field fluctuations yields 
an equation that is very similar to equation (Q) except 
that our function f{K) — {fi) due to phase fluctuations 
is replaced by the expectation value {A"^). For d = 3, 
magnetic fluctuations of the gauge field A produce a first 
oder transition and move the critical point into the posi- 
tive cr region. The domain of validity is also in the small 
cr region as for our approach. 

We now show that the origin of the first order tran- 
sition of equations ( pT| ) and (|l^) is not due to the ap- 
proximation itself and that the domain of validity of the 
approximation is precisely in the first order domain for 
d = 3. We establish a quantitative criterion by com- 
paring the size of the critical gaussian region according 
to the Ginzburg criterion with the critical temperature 
due to phase fluctuations. The Ginzburg criterion |l^] 
measures the importance of correlated gaussian fluctua- 
tions that are mainly due to the amplitude, i.e. they are 
negligible if 



\t\ » ta 



2k 

2Vo 



where tc deflnes the critical region. \tc\ from equation (^ 
is a measure of the importance of phase fluctuations. Our 
approach is then valid when (amplitude) fluctuations are 
dominated by phase fluctuations, i.e. when \tc\ >> Iq. 
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Setting uj ^ d—1 which is almost the correct XY critical 
temperature, we get the criterion: 

(13) 

For d = 1, the criterion is in agreement with our results 
(not shown here) which are almost identical to the exact 
solution jl^. For d — 2, the criterion is not clear. The 
approximation is reliable precisely in the first order do- 
main for the 3d case (uTc/Vo << 1), and everywhere for 
d = 4. For d > A, the domain of validity is in the contin- 
uous domain. Therefore, we conclude that the first order 
transition is not a consequence of the approximation for 
d = 3, whereas a doubt remains for d = 2. 

The first order transition is due to the fact that the en- 
ergy of vortices can be lowered by a reduced value of the 
amplitude, which is energetically favorable when the po- 
tential energy of the amplitude is sufficiently soft. In this 
case, the usual transition scenario, given by an unbind- 
ing of topological excitations, seems to be replaced by a 
sudden proliferation of the latter. Such a process may 
be difficult to describe within the renormalisation group 
e-expansion where the interplay between amplitude and 
phase is not explicitely taken into account [ pT| . A similar 
scenario has been found by Minnhagen et al. ^ for the 
2d Coulomb gas. They show that a large vortex fugacity 
y — exp(— /37|'0p7r^/2) can lead to a discontinuous tran- 
sition produced by a proliferation of vortices, whereas 
small fugacity causes the usual Kostcrlitz-Thouless tran- 
sition. When (T is small, our scenario corresponds to the 
case of large vortex fugacity. The mean amplitude can 
become much smaller than the one of the XY model (see 
figure 1 where R~l corresponds to the XY model) . For 
large ct (i? w 1) , the KT transition is recovered. Our 
result (ctc ~ 1.25) is in good agreement with the one of 
Bormann and Beck 1^ who found a first order transition 
for CT < 1. 

For real superconductors, a is smaller than 1 
Therefore, superconductors with low Tc/Vq could have 
a first order transition. A candidate for a possible ob- 
servation would be a 3d superconductor with a critical 
region to that is not too small but still in the first oder 
domain. In the BCS limit, the transition becomes very 
weakly first order such that there is no contradiction with 
the continuous behaviour of BCS superconductors. Un- 
derdoped cuprates are quasi 2 dimensional and have a 
low Vb/Tc. Their transition is then XY-like, whereas 
overdoped cuprates, that are almost 3d and have a large 
Vq/Tc, could have an observable first order transition. 
It is, however, interesting to remark that measurements 
of the entropy change AS* at the vortex lattice melting 
transition of Bi2Sr2CaCu208 |jl^ have shown a dramatic 
increase in AS* per vortex when the zero field transition 
is approached. This could be a hint that the supercon- 
ducting transition remains first order even in zero field. 



In this paper we have investigated the thermodynamic 
properties of the classical Ginzburg-Landau (GL) model. 
It is determined by two model parameters, tr and Vq. 
a governs the strength of amplitude fluctuations and Vq 
the overall strength of fluctuations of the complex GL 
field. We have treated the model by a variational ap- 
proximation which takes into account the coupling be- 
tween phase and amplitude through effective coupling 
constants. Minimizing the corresponding variational free 
energy leads to a set of self-consistent modified GL equa- 
tions containing phase and amplitude fluctuations. The 
behaviour of the GL transition changes when the ratio 
cr/(Vo/Tc) is varied: for cr < C Vq/Tc the transition of 
first order with C = 9.9 for d = 3 and C = 2.2 for d = 2. 

For c? = 3, we showed that phase fluctuations dominate 
the transition in the first order domain and amplitude 
fluctuations can be neglected. We thus conclude that 
a flrst order transition is indeed a valid scenario for the 
GL model, once the amplitude of the fleld can sufficiently 
adapt in order to lower the total energy of the system. 
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